Physical origin of the power-law tailed statistical distributions 
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Starting from the BBGKY hierarchy, describing the kinetics of nonUnear particle system, we 
obtain the relevant entropy and stationary distribution function. Subsequently, by employing the 
Lorentz transformations we propose the relativistic generalization of the exponential and logarithmic 
functions. The related particle distribution and entropy represents the relativistic extension of the 
classical Maxwell-Boltzmann distribution and of the Boltzmann entropy respectively and define 
the statistical mechanics presented in [Phys. Rev. E 66, 056125 (2002)] and [Phys. Rev. E 72, 
036108 (2005). The achievements of the present effort, support the idea that the experimentally 
observed power law tailed statistical distributions in plasma physics, are enforced by the relativistic 
microscopic particle dynamics. 
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I. INTRODUCTION 

In plasma physics, the power-law tails in the particle 
population, has been systematically observed in the last 
fifty years. For instance the cosmic ray spectrum 



(1) 



which extends over 13 decades in energy, from 10^ eV 
to 10^° eV, and spans 33 decades in particle flux, from 
10~^^ to 10^ units, obeys the Boltzmann law of classical 
statistical mechanics i.e. 



X{x) exp(-x) 



(2) 



for low energies, while for high energies this spectrum 
presents power law fat tails i.e 



X{x) - X' 



(3) 



the spectral index k being close to 0.32-0.37. 

The above spectrum was approached for the first time 
in 1968, by using a different distribution from the Boltz- 
mann one. In his proposal Vasyliunas heuristically iden- 
tified the function xi^) with the Student distribution 
function which presents power-law tails [l|. In the last 
40 years a vast amount of literature has been produced, 
regarding the so called kappa-plasmas, based on the Va- 
syliunas distribution. Up to now, several physical mech- 
anisms have been explored in order to furnish theoretical 
support to the experimentally observed power-law-tailed 
distribution functions. However there is currently an in- 
tense debate regarding the theoretical foundations of the 
non-Boltzmannian distributions. 

In the last years, after noting that the power-law tails 
are placed in the high energy region, and then regards rel- 
ativistic particles, the question has been posed whether 
the solution of the problem, i.e. the theoretic determina- 
tion of the function xi^) ^^'^ consequently of the related 
distribution and entropy, can be explained by invoking 
the basic principles of special relativity. 



The present paper, deals with the relativistic statistical 
theory [2l-[l0j. predicting for the function x{^): the very 
simple form i.e. 



X{x) = exp^(-x) 



with 



exp^(x) = (\/ 1 + K^x'^ -f KxY^'^ 



(4) 



(5) 



The parameter k < 1 is the reciprocal of light speed in a 
dimensionless form while the function exp^ (x) represents 
the relativistic generalization of the ordinary exponential 
which recovers in the classical limit k — > 0. 

In the last few years various authors have considered 
the foundations of the statistical theory based on the dis- 
tribution function involving the generalized exponential 
exp^(a;), e.g. the H-theorem and the molecular chaos 
hypothesis [ill . [T2j , the thermodynamic stability [l^ [ij] , 
the Lesche stability ^,5- 18J, the Legendre structure of the 
ensued thermodynamics Il9i. 20l | , the thermodynamics of 
non-equilibrium systems [21 1, quantum versions of the 
theory [22|-|2^. the geometrical structure of the th eory 
p5l |. various mathematical aspects of the theory [26l - l3ll |. 
etc. On the other hand specific applications to physi- 
cal systems have been considered, e.g. the cosmic rays 
y , relativistic [32] and classical fSS"] plasmas in presence 
of external electromagnetic fields, the relaxation in rela- 
tivistic plasmas under wave-particle interactions 
anomalous diffusion [sgI [37j , non-linear kinetics 
kinetics of interacting atoms and photons [13], particle 
kinetics in the presence of temperature gradients 
particle systems in external conservative force fields 
stellar distributions in astrophysics [isl - lisj . quark-gluon 
plasma formation (46| , quantum hadrodynamics models 
|47j . the fracture propagation [i^, etc. Other applica- 
tions regard dynamical systems at the edge of chaos [4§- 
'51'], fractal systems (s^, field theories fss f, t he random 
matrix theory js^ . l55l]. the error theorv|56| , the game 
theory [s^l, the information theory [sl], etc. Also appli- 
cations to economic systems have been considered e.g. to 
study the personal income distribution to model 

deterministic heterogeneity in tastes and product differ- 
entiation |63| etc. Finally in [6^1, some historical 
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remarks, on the theories deahng with power-law tailed 
distribution functions, are reported. 

In the present contribution we reconsider critically the 
foundations of the statistical theory the generalized ex- 
ponential (O. Our main goal is to show that (i) the 
function exp^ (x) , can be obtained within the one-particle 
relativistic dynamics, in a very simple and transparent 
way, starting directly from the Lorentz transformations 
(ii) the standard principles of ordinary relativistic kinet- 
ics, conduct unambiguously to the relativistic general- 
ization of the classical Boltzmann entropy and Maxwell- 
Boltzmann distribution, without invoking any extra prin- 
ciple or assumption. 



II. KINETIC EQUATION 

Let us consider the most general relativistic equation 
imposing the particle conservation during the evolution 
of a many body system. That equation is the first equa- 
tion of the BBGKY hierarchy i.e. 



p^d^f-mF 



dp"" 



P'° Pi P'l 



d^p' (fpi d^p[ ^ „, „ . , 

(6) 

and describes, through the one-particle correlation func- 
tion or distribution function / — f{x,p), a relativistic 
particle system in presence of an external force field. The 
streaming term as well as the Lorentz invariant integra- 
tions in the collision integral, have the standard forms of 
the relativistic kinetic theory [6^, [g^I . 

The two particle correlation function [6^ , here denoted 
by f®fi, appearing in Eq. ((5]), at the moment remains 
an unknown, two variable function. We recall that in 
classical kinetics, the two-particle correlation function, 
according to the molecular chaos hypothesis, is the ordi- 
nary product of / and /i i.e. f ® fi = f fi- Therefore 
the composition law f ® fi can be viewed as a relativis- 
tic generalized product of / and /i, isomorphic to the 
ordinary product. 

Following standard lines of kinetic theory, we note that 
in stationary conditions, the collision integral in Eq. Q 
vanishes and then it follows that f®fi — f ®fi- More 
in general it holds 



L(/®/i) = L{f'®f[) , 



(7) 



L{x) being an arbitrary function, and this relationship 
expresses a conservation law for the particle system. On 
the other hand a conservation law has form 



A(/) + A(/i) = A(/') + A(/0 



(8) 



A(/) being the collision invariant of the system. There- 
fore we can pose 



L(/®/i) = A(/) + A(/i) . 



(9) 



From the definition of the correlation function and tak- 
ing into account of the indistinguishability of the parti- 
cles, it follows that f®\ = 1(E) f = / and this implies that 



A(l) = and L{f) = A(/). After taking into account 
the later relationship, Eq. © becomes 



A(/®/i) -A(/) + A(/i) 



(10) 



Consequently the collision invariant A(/) permits us to 
determine univocally the correlation function f ® fi as 
follows 



/®/i=A-i(a(/)+A(/i)) 



(11) 



In relativistic kinetics, the collision invariant A(f), un- 
less an additive constant, is proportional to the micro- 
scopic relativistic invariant /. Then we can pose 



A{f) = -^I + I3fi 



(12) 



(3 and /3/i being two arbitrary constants. In presence of 
an external electromagnetic field A'^ , the more general 
microscopic relativistic invariant /, has a form propor- 
tional to 



/ = {p" + qA^/c) - mc^ , 



(13) 



C/y being the hydrodynamic four-vector velocity with 

VU^ = Finally, after inversion of Eq. p^ . 
the stationary distribution is obtained as follows 

f = K-\-PI + P^x) . (14) 



It is remarkable to note that Eq. (|12p follows from the 
variational equation 



5_ 

Sf 



S 



d^plf + f3fi d^pf 



, 



(15) 



where the functional 5, unless an arbitrary additive con- 
stant, is given by 



S 



d^p 



A(/) df 



(16) 



/ A(/) df being the indefinite integral of A(/). The vari- 
ational equation (1151) represents the maximum entropy 
principle. The constants (3 and l3fj, are the Lagrange mul- 
tipliers while the functional S, defined though Eq. p6|) . 
is the system entropy. 

We stress that the function A(/) defines univocally 
both the stationary distribution ([T^ and the entropy 
p6| of the system as well as the two-particle correlation 
function pT|) . In classical statistical mechanics, it is well 
known that A(/) = ln(/) so that the two particle correla- 
tion function becomes / <8) /i = / /i , while and (fTB]) 
reduces to the exponential distribution and Boltzmann 
entropy respectively. 

In the next section, in order to develop a relativistic 
statistical mechanics, we will determine the function A(/) 
within the special relativity, starting from the Lorentz 
transformations. 
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III. LORENTZ TRANSFORMATIONS 

In the present section we will show that the function 
^~^{f) emerges within the special relativity as the rela- 
tivistic generalization of the ordinary exponential of clas- 
sical physics. 

Let us consider in the one-dimension frame S two iden- 
tical particles A and B, of rest mass m. We suppose that 
the two particles move with velocity Va and Vb respec- 
tively. The momenta of the two particles are indicated 
with Pa = p{va) and Pb = p{vb), while their energies 
are indicated with Ea — E (va) and Eb = E (vb) respec- 
tively. 

In classical physics, in the rest frame S' of particle 
B, the momentum of the particle B is p^ — while 
the momentum of the particle A is given by the Galilei 
transformation formula 



Pa 



Pa 



Pb 



(17) 



After introducing in place of p the dimensionless mo- 
mentum q = pIp*, we note that the exponential function 
exp((7), permits us to write the Galilei additivity law (1171) . 
in the following factorized form 



exp(g^) = exp(g^)exp(-(7B) 



(18) 



The Galilei relativity principle, imposes the equiva- 
lence of all the inertial frames. According to this prin- 
ciple, the inverse Galilei transformation must have the 
same structure of the direct transformation ([T51) except 
for the substitutions O and <7b — >■ — gs- This re- 
quirement is satisfied thanks to the following property of 
the exponential function 



exp(a;) exp(— x) = 1 



(19) 



We consider now the two particles in the rest frame S' 
of particle -B, within the special relativity. The velocity, 
momentum and energy of the particle B are = 0, 
p'b = and E'^ = mc^ respectively. In S' the velocity of 
particle A is given by the formula 



(20) 



1 - VaVb/c^ 



defining the relativistic velocity composition law, which 
follows directly from the kinematic Lorentz transforma- 
tions. In the same frame S' the momentum and energy 
of particle A are given by the dynamic Lorentz transfor- 
mations 

Pa = 1{Vb)Pa - C~'^VbJ{Vb)Ea , (21) 

E'a = j{vb)Ea - Wb7("b)?'a , (22) 

i{vb) = (1 — ^b/c^)~^^^ being the Lorentz factor. Af- 
ter taking into account the expression of the momentum 
Pb — mvB^ivB) and of the energy Eb — mc^^{vB) of the 
particle B the latter transformations become 

Pa = PAEB/mc^ - EAPBlmt? , (23) 
E'^ = EaEbItu^ - PaPbItu , (24) 



Let us introduce in place of the dimensional variables 
{v,p,E) the dimensionless variables {u,q,£) through 



V 

u 




|k|c: 



< c 



(25) 



From its definition, it follows that k can be viewed 
as the Einstein /3 factor related to the velocity w*. The 
condition < c, implies that — 1 < k < +1. 

The dynamic Lorentz transformations for the dimen- 
sionless momentum and energy variables become 



^A — ^ Qa£b Qb£a I 

£a " ^ £a£b ^ QaQb ; 



(26) 
(27) 



while the classical limit c — > oo is replaced by the limit 

For a particle at rest it results E{0) — mc^ and then 
£"(0) = 1/k^. Then represents the dimensionless 

rest energy of the particle. Alternatively 1/k can be 
viewed as the refractive index of a medium in which the 
light speed is w*. 

From the Lorentz invariance it follows easily the 
energy-momentum dispersion relation 



4c2 



2 2 1 

K q =1 



(28) 



After expressing in the right hand side of Eq. ((26l) . the 
energy in terms of the momentum £ = -y/l+Tt^g^/K^, 
we obtain the momentum relativistic additivity law as 
follows 



q 



A = qA^/TTK^ - qB\/T+l^ql , 



(29) 



which in the classical limit /t — )• reproduces the ordinary 
additivity law ([T7|) . of classical physics. 

The Galilei relativity principle, holding both in classi- 
cal physics and in special relativity, imposes the equiva- 
lence of all the inertial frames. According to this prin- 
ciple, the inverse transformation of (|29|) must have the 
same structure of the direct transformations (1^^ except 
for the substitutions O Qa and g^ — ^ —qs i.e. 



Qa = + i^'^ll +1b\/i + k'^Qa'^ , 



(30) 



It is easy to verify that ([501 follows directly from ([^^ 
and viceversa. 

The Lorentz transformation for the relativistic mo- 
menta (|29|) . representing the momenta additivity law in 
special relativity, has the important feature that the con- 
tributions of the two particles, appearing in the right 
hand side of the equation, are not factorized. Sponta- 
neously the question emerges at this point, whether new 
variables exist, able to factorize the contribution of the 
two particles in the right hand side of the relativistic ad- 
ditivity law ([29]). 

It is easy to verify that the new variable is given by 
function 



exp«(g) = ( VI + kV 



Kg 



1/k 



(31) 
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so that the Lorentz transformation 
fowing factorized form 



29l) assumes the fol- 



(32) 



with 



On the other hand Gahlei relativity principle imposes, 
for the inverse Lorentz transformation, the following fac- 
torized form 



exp«(gA) = exp^(gB) exp^{q[J 



(33) 



By comparison of the above direct and inverse Lorentz 
transformations, it obtains the property 



exp«(9)exp^(-(7) = 1 , 



(34) 



which can be verified easily, by direct inspection of the 
definition ([3T|) . 

It is remarkable that the function exp^(q) emerges in 
one-particle special relativity as the variable, able to fac- 
torize the momentum Lorentz transformation, and repre- 
sents the relativistic generalization of the ordinary expo- 
nential which factorize the momentum Galilei transfor- 
mation in classical physics. Clearly in the classical limit 
K — >■ 0, exp^(5) reduces to exp(g). 

The inverse function of exp^ (g) indicated by ln„ (q) and 
defined through ln„(exp^((7)) = exp^(lnK((?)) = 1, repre- 
sents the relativistic generalization of the ordinary loga- 
rithm, which recovers in the classical limit k — > 0, and is 
given by 



2k 



(35) 



The property of exp^(q), enforced by the Galilei rel- 
ativity principle, transforms into the following property 
of In^ix) 



(36) 



holding also for the ordinary logarithm of classical 
physics. 



IV. RELATIVISTIC STATISTICAL 
MECHANICS 

In the previous section it has been shown that the func- 
tions exp^(g) and InK(g) emerge as the relativistic gen- 
eralizations of the ordinary exponential and logarithm 
functions of classical physics. Therefore in the following 
we pose 



A(x) = lnf,{x) , 
A"i(a;) = cxp„(a;) . 

After taking into account of the property 

X ln„(a;) = - ln^,{ex) , 
ax 7 



(37) 
(38) 



(39) 



1 



7 = 



vr 



K 



2 ' 



e = exp^(7), 
the entropy (fT6)) simplifies as 

j3„ 



S 



(40) 
(41) 

(42) 



It is worth stressing that the latter relationship de- 
fines the relativistic entropy, as proportional to the mean 
value of — lnK(//e), like in the case of classical sta- 
tistical mechanics where the Boltzmann entropy, S — 

— J (Pp f In (//e), is proportional to the mean value of 

— ln(//e). Clearly the entropy (H^ . in the classical limit 
reduces to the Boltzmann entropy. The constant 7, given 
by represents the Lorentz factor related to the ve- 
locity appearing in (j25p . and in the classical limit, 
approaches the unity. On the other hand the constant 
e given by (j41l) . represent a relativistic generalization of 
the Napier number e, which recovers in the classical limit. 

The entropy (UU can be written explicitly as follows 



2k 



+ K. 



1 



1 



(43) 



while the related stationary distribution (|T4|) assumes the 
form 



/ = exp„(-/3/ + /?/i; 



(44) 



and reduces to the Maxwell-Boltzmann distribution, in 
the classical k — limit. 

The distribution (|44p . in the global rest frame where 
Ui, = (c, 0, 0, 0) and in absence of external forces i.e. 
A'^ — 0, simplifies as 



/ = exp, (-PE + p^^) 



(45) 



E being the relativistic kinetic energy. This distribution 
at low energies {E — )■ 0) reduces to the classical Maxwell- 
Boltzmann distribution i.e. / oc exp(— /3i?), while at 
relativistic energies {E — > +00) presents power-law tails 
i.e. 



f^E 



(46) 



in accordance with the experimental evidence. We recall 
that the first experimental validation of the distribution 
()45p . concerns cosmic rays and has been considered in 
ref. [1]. Recently a computer validation of the same 
distribution, has been considered in refs. [13, HI], where 
the relaxation in relativistic plasmas under wave-particle 
interaction, has been simulated numerically. 

Finally, after posing A(/) = lnB;(/), in Eq. ((TT|) . the 
relativistic two-particle correlation function assumes the 
form 

=exp,(ln,/ + ln«/i) , (47) 

which in the k — >■ classical limit, reduces to //i, as dic- 
tated by the molecular chaos hypothesis. Consequently 
the relationship (|T7)) can be viewed as defining the rela- 
tivistic extension of the molecular chaos hypothesis. 
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V. CONCLUSIONS 

Let us consider the probability distribution function 

/ = eiexp,(-6/3[£;-M]) , (48) 

^1 and ^2 being two arbitrary constants. In the expres- 
sion of the distribution function E is the relativistic mi- 
croscopic energy, while /? and {3^ are the Lagrange mul- 
tipliers. The Maximum Entropy Principle asserts that 
the distribution (j48p can be obtained, by maximizing the 
entropy 

S=-^J d'p flnM/e^i) , (49) 

with 7 — l/Vl — and e = exp^(7), under the con- 
straints imposing the conservation of the norm of / and 
of the mean value of the E. It is remarkable the above 
distribution and entropy are linked through the Maxi- 



mum Entropy Principle, independently on the particular 
values of the arbitrary constants and ^2 • 
In kinetic theory, customarily it is posed 

a = 6 = 1 , (50) 

in order to simplify the expression of the distribution 
function and of the two particle correlation function, ap- 
pearing in kinetic equation. This choice, is made natu- 
rally, in the present paper, being the starting point of our 
presentation the kinetic equation ^ . 

On the other hand, in statistical mechanics, in order 
to simplify the expression of the entropy, it is posed 

= 1/e , (51) 
6=7- (52) 

This latter choice has been made for instance in the ref. 
[1, [1] where the starting point for the presentation of the 
theory was the entropy functional. 
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